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ABSTRACT 

Context. The pioneer space mission for photometric planet searches, CoRoT, steadily monitors about 12,000 stars in each of its 
fields of view; it is able to detect transit candidates early in the processing of the data and before the end of a run. 
Aims. We report the detection of the first planet discovered by CoRoT and characterizing it with the help of follow-up observations. 
Methods. Raw data were filtered from outliers and residuals at the orbital period of the satellite. The orbital parameters and the 
radius of the planet were estimated by best fitting the phase folded light curve with 34 successive transits. Doppler measurements 
with the SOPHIE spectrograph permitted us to secure the detection and to estimate the planet mass. ** 

Results. The accuracy of the data is very high with a dispersion in the 2.17 min binned phase-folded light curve that does not 
exceed ~ 3. x 10~ 4 in flux unit. The planet orbits a mildly metal-poor GOV star of magnitude V=13.6 in 1.5 days. The estimated 
mass and radius of the star are O.95±O.15M0 and 1.11±O.O57?0 . We find the planet has a radius of 1.49±0.087?.j up , a mass of 1.03 
±0.12Mj up , and a particularly low mean density of 0.38 ±0.05g cm -3 . 

Key words, planetary systems - techniques: photometry - techniques: radial velocity 



1. Introduction 

CoRoT is the first space survey dedicated to the photomet- 
ric search for extrasolar planets. Led by the CNES space 
agency, CoRoT was born from a joi nt effort of France , 
European countries, ESA, and Brazil (jBaglin et al. 2 006b 
The detailed description of the instru ment and the mission 
can be found in a pre-launch book (|CoRoT 2006f ) and in 
a post-launch paper (Auvergne et al. in preparation). The 
instrument satisfies the initial scientific requirements of the 
mission with a noise level that matches the photon noise 
over most of the magnitude range accessible to CoRoT (12- 
16). The measure of the stellar flux every 512 s results from 
the piling up on board of 16 individual exposures of 32 s 
obtained by aperture photometry. The photometric masks 
on the CCDs are uploaded at the beginning of a run, af- 
ter being selected from a pre-defined list of patterns. In 
addition to the target, a mask may contain contaminating 
stars that are potential sources of noise and false alarms. 
Reference windows with a square shape were also selected in 
the non-exposed parts of the CCDs (5x5 pixels) the offset 
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* The CoRoT space mission, launched on Dec. 27th, 2006, was 
developed and is operated by the CNES, with participation of 
the Science Program of ESA, ESTEC/RSSD, Austria, Belgium, 
Brazil, Germany and Spain. Based in part on observations with 
the SOPHIE spectrograph at Obs. de Haute Provence, France 

** Table[TJis available in electronic form at the CDS via anony- 
mous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or the web site. 



corrections or in the darkest regions of the CCDs (10 x 10 
pixels) for sky-background corrections. In polar Earth or- 
bit, CoRoT is free of the main limitations of ground-based 
transit searches. The photometric accuracy is better than 
lmmag and the monitoring of the targets is nearly contin- 
uous over several months. However, the total duration of a 
run cannot exceed 6 months to avoid blinding by the Sun 
and occultation by the Earth. Another inconvenience comes 
from charged particle impacts at the crossing of the South- 
Atlantic Anomaly (SAA) as they can produce bright pixels 
on the CCDs and outliers in the data. This letter reports 
the detection and characterization of the first planet dis- 
covered by CoRoT in its initial run of observation. With a 
large radius and a low mean density, the planet seems to de- 
part from the distribution of the known transiting planets. 
The results are based on a preliminary analysis of the raw 
data that we briefly describe and on ground-based follow- 
up with the SOPHIE spectrograph. The full analysis of the 
light curve will be presented elsewhere using the fully re- 
duced data issued from the standard pipeline of corrections. 

2. CoRoT observations 

CoRoT's observation sequence started with a pointing di- 
rection close to the anticenter of the Galaxy. This initial run 
lasted from February 6th to April 2nd 2007 and supplied 
~12,000 high accuracy light curves for exoplanet transit 
search. These light curves are nearly continuous over 55 
days with only a small number of gaps that mainly result 
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from the crossing of the SAA (the duty cycle is 92%). The 
nominal sampling is 512 s but can be 32 s for —10% of 
the targets (500 out of 5700 in each of the two CCDs). A 
peculiarity of the CoRoT exoplanet program is the "alarm 
mode" , an operational loop between the science team and 
"Control & Command Center" on ground and the instru- 
ment in space. Its goal is to optimize the science return 
of the mission by identifying transit candidates early in 
the process, before the end of run and before the data are 
fully reduced. For these transit candidates a decision can be 
made to change the rate of the observations from one expo- 
sure in 512 s to one exposure in 32 s. The interest for this 
oversampling is two-fold: (i) to get a better coverage of the 
transit profile; (ii) to reduce the level of noise by removing 
outliers from the corrupted exposures. Last but not least, 
this "on the stream" detection makes it possible to start 
follow-up operations as soon as possible. Thanks to the 
alarm mode, some transit candidates were identified soon 
in the raw data and the sampling rate changed correspond- 
ingly. For example, the best candidate, CoRoT-Exo-lb, 
found around a GOV star of 7=13.6 (CoRoT-Exo-1) with 
a period of 1.5 days and a depth of 2%, was submitted to 
a series of follow-up observations that started at mid run. 

3. Data and preliminary treatments 

The light curve of CoRoT-Exo-1 consists of -68,824 data 
points for a total duration of 54.72 days. Sampling is 512 s 
during the first 30 days and 32 s for the rest of the time. 
Out of 36 detected transits, 20 are found in the part of 
the light curve with nominal sampling and 16 are found 
in the oversampled part. Data used for "alarm" detection 
are raw data requiring further processing to get correct es- 
timates of the transit parameters, particularly, the transit 
depth. We used a number of treatments that focus: (i) on 
the background and offset corrections; (ii) on the filtering of 
the outliers and low frequency residuals produced at SAA 
crossings. The developed procedures are forerunners of the 
ones implemented in the reduction pipeline. 
— Reference windows for the background and the offset 
corrections can be affected by charged particle impacts at 
SAA crossing. To detect and eliminate, in these windows, 
the resulting outliers, we subtracted a moving-median fil- 
tered version of the light-curve and discarded the points 
at distances greater than a few times the dispersion of the 
residuals. For the background, a constant level was esti- 
mated using the reference sky window closest to the tar- 
get. The background light curve, once filtered from the 
outliers, is smoothed out with a 190 points width me- 
dian filter and interpolated into the data points (as the 
time sampling changed during the observations). By this 
method the average background level is 11.9 e"px~ 1 s _1 
(5.58 ADU px~ 1 s _1 ). For the offset, we proceeded in the 
same way as for the background. The camera offset is ob- 
tained as the median of 23 reference windows. The offset 
curve is low-pass filtered with the same median filter as 
above, but with a width of 500 points. Thus, the median 
value of the offset is 55.47e^px- 1 s" 1 (25.9 ADU px^V 1 ). 
To get the total correction, we made the sum of the back- 
ground and the offset (estimated within the photomet- 
ric mask) and subtracted the result to the original data. 
After multiplying by the detector gain, we found the me- 
dian flux from the star is 221,000e~ per 32 s exposure 
(640.45e _ px _1 s _1 ) with a signal-to-noise ratio of 476.19 . 



— The CoRoT-Exo-1 light curve also contains outliers that 
are identified and located thanks to a median filtering, in 
the same way as performed for the background and off- 
set light curves. With such a procedure, 8% of the data 
points were identified as outliers and eliminated from the 
subsequent analysis. The filtered light curve keeps the sig- 
nature of the satellite orbital period (at 102 min), likely 
due to residuals of the background corrections. We modeled 
this effect by estimating, for each orbit of the satellite, the 
modulation obtained from the 50 closest orbits; the light 
curve is successively: (i) normalized with the median value 
at each orbit, (ii) phase- folded at the period of the satel- 
lite, (iii) high frequency filtered (using a 128 x 128 points 
Savitzky-Golay weighted moving average filter, with 4th 
degree weighting). From the preparatory ExoDat database 
only two faint contaminants (V > 17) are present within 
the 61 pixels of CoRoT-Exo-l's photometric mask; their 
contribution to the total flux was neglected. 

4. Folded light curve and best fitting 

Possibly due to bright pixels inside the photometric mask, 
the light curve has low frequency (< lcycle / day) residuals 
that could mar the transit shape and affect the estimate 
of the star and planet parameters. To minimize such effect 
we proceeded as follows: (i) around each individual tran- 
sit (just before and after), a parabola is chosen to fit the 
data points close to the transit phase, (ii) the fit is high 
frequency filtered and subtracted to the ensemble of points 
ranging from the transit orbital phase —0.1 to the phase 0.1, 
(iii) the 36 transits are re-sorted in orbital phase, using the 
ephemeris reported in Tabled Then, the points are binned 
in 0.001 phase- wide bins (2.17 min in time unit), and the er- 
ror bars are estimated as the dispersion of the points inside 
the bin divided by the square root of the number of points 
per bin. The median number of points per bin is 316 and the 
average la error bar estimated this way is 0.0003 (~ 3 times 
larger than expected from the photon noise limit). The 
noise level should be reduced once the pipeline includes jit- 
ter corrections. The detrended light curve was phase- folded, 
avoiding 2 transits strongly affected by a bright pixel (see 
Figure [T]), and fitted to a model following the formalism of 
Gimenez (2006). To find the solution that best matches the 
data, we minimize d the \ 2 using the algorithm AMOEBA 
([Press et al. 19 92). The fitted parameters are the transit 
center, the orbital phase at transit start, Rp/R+, the or- 
bital inclination i and the two non linear limb-darkening 
coefficients u + and u_ (jBrown et al. 200~iT ). The lcr uncer- 
tainties were estimated by a bootstrap analysis based on 
simulated data sets constructed from the best fit residuals. 
Building each data set required: (i) to subtract the best fit 
model to the data, (ii) to sort again a fraction of 1/e of the 
residuals, (iii) to add again the best fit model to the new 
residuals. AMOEBA was used to fit all the parameters in 
each of the new data sets. The uncertainties were estimated 
as the standard deviation of each of the fitted parameters 
(see Table [2). 

5. Follow-up observations 

Observations of CoRoT-Exo-1 were reported in the archive 
of the pre-launch photometric survey performed with the 
Berlin Exoplanet Search Telescope (BEST); the data how- 
ever, were too noisy to clearly identify a transit egress. 
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Fig. 2. Radial-velocity variations of CoRoT-Exo-lb versus 
phase from CoRoT's ephemeris. Top: the data fitted with 
a circular orbit of semi-amplitude K=188ms~ 1 and a drift 
of 1.02 ms -1 per day; bottom: the O-C residuals to the fit. 



Fig. 1. Normalized and phase folded light curve of the 34 
best transits of CoRoT-Exo-lb (top); the residuals from 
the best- fit model (bottom). The bin size corresponds to 
2.17mm, and the la error bars have been estimated from 
the dispersion of the data points inside each bin. 



Follow-up photometry with the Wise observatory 1.0 m 
telescope confirmed that a transit occurs on the main 
target. Images taken at the CFHT showed only some 
faint background stars (up to V=22) that confirm the 
weak contamination inside CoRoT's photometric mask. 
The results of these investigations, indicative that none 
of the CoRoT-Exo-l's near neighbour can be responsi- 
ble for the transit signal, will be published elsewhere. 
High precision radial velocity observations of CoRoT-Exo-1 
were made at Observatoire de Haute Provence in March- 
April and October 2007 with SOPHIE, an echelle cross- 
dispersed fiber-fed spectrograph at the 193 cm telescope 
(|Bouchv et al. 20061 ). The spectrograph routinely gets to a 
precisi on less than 2-3 ms" 1 with the cross-correlation tech- 
nique (|Baranne et al. 1 996) on bright targets. CoRoT-Exo- 
1 was observed with the High Efficiency mode (R=40,000) 
at an accuracy of about 20-30 ms -1 in typical exposure 
times of 45 min. Data reduction was made online and cross- 
correlation performed with a mask corresponding to a G2 
star. The observations polluted by Moon illumination have 
been corrected for the correlation peak of the sky back- 
ground, which is measured simultaneously through a neigh- 
bour fiber. The estimated error bars take into account such 
residual systematics that affect the Doppler measurements 
on faint targets, with signal-to-noise ratios ranging from 
20 to 27. In total, 9 measurements were performed with 
SOPHIE and are reported in Table Q] where are displayed: 
Julian date, planetary orbital phase, radial velocity and er- 
ror. The measurements show a variation in phase with the 
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Fig. 3. Bisector spans: no variation with amplitude 
and phase corre lates with the radial velocities (e.g. 
lOueloz et al 200l[) 



ephemeris constrained, with a very high accuracy, by the 
CoRoT light curve. The orbital period and the epoch of 
the transits are fixed to the CoRoT values. As the period is 
very short, tidal circularization justifies a zero eccentricity 
assumption. The semi-amplitude of the radial velocity vari- 
ation, K =188±llms _1 , is compatible with a companion 
of planetary mass and was adjusted to the data (see Figure 
[2]). The best fit is obtained assuming a drift of 200 ms -1 be- 
tween the two epochs of SOPHIE's observations or, equiva- 
lently, a linear drift of 1 ms -1 per day. The final solution is 
displayed on Tabled the O-C residuals have a standard de- 
viation of 34 ms -1 . The measurements do not correlate with 
the bisector of the cross-correlation function (see Figured]). 
In summary, the radial-velocity curve supports the plane- 
tary nature of the transiting body detected by CoRoT and 
discards other interpretations involving background stars, 
grazing eclipsing binary or a triple system. Additional ob- 
servations are ongoing with HARPS to clarify the nature 
of the drift, possibly due to a second companion. 
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6. Star and planet parameters 

To derive the star's parameters, we took advantage of the 
spectroscopic observations performed to get the radial ve- 
locity measurements. The very first observations recorded 
with HARPS, whose resolving power is R ~ 100,000, 
were used to perform a spectral analysis. For each set 
of spectra, we co- added the individual exposures, order 
per order, once calibrated in absolute wavelength and re- 
binned at a constant step of 0.02A. The resulting signal- 
to-noise ratio per pixel at 550 nm is ~ 80 for the co- 
added HARPS spectra. The atmospheric parameters, T g, 
log g and [M/H] were derived by comparing the observed 
normalized spectrum to a library of synthetic spectra ca l- 
culated with LTE MARCS models (jGustafsson et al. 20051) . 
The library sampled the range of temperature from 4000 K 
and 8000 K, surface gravity from log g = 1.0 to 5.0 and 
metallicity from —3.0 up to 5.0. The temperature was first 
estimated on the Ha line wings, which is insensitive to any 
stellar parameters but the temperature. Special care has 
been devoted to the normalization of the corresponding or- 
der which could be quite tricky due to the extended wings 
of line. For this line, we u sed the new grids of ATLAS9 
models ()Heiter et al. 20 02). This value was controlled us- 
ing the relative intensity of iron lines, from which a simi- 
lar temperature was derived. Besides, this is also in good 
agreement with that estimated from its color index in the 
ExoDat database (B-V = 0.573) and the standard photo- 
metric calibrations. Metallicity is difficult to determine in 
this part of the spectrum as all indicators are also sensi- 
tive to gravity. Both were estimated jointly by comparing 
the observed spectrum with the synthetic ones, until the 
best fit was achieved. In our modeling of the spectra, we 
used the rotational velocity vaini =5.2±1.0kms _1 derived 
from the widths of th e HARPS cross-correlation func- 
tion (jSantos et al. 2002] ). Though more detailed studies are 
needed to achieve a precise determination, we find that 
CoRoT-Exo-1 is similar to the Sun with a temperature of 
T cff = 5950K ±150 K, log g =4.25±0.30 but a metal defi- 
ciency [M/H]=— 0.30 ± 0.25. This metallicity is in agree- 
ment with the [Fe/H] value of —0.4 deduced from the 
calibration of the H ARPS cross-correla tion function and 
the star's color index (jSantos et al. 2002f ). The stellar mass 
was derived using the evolutionary tracks (jSiess et al. 2 000; 

iMorel fc Lebreton 20071) in the (Ml /3 /R*,T eS ) plane. The 

1 /3 

method allows us to take advantage of the M* / i?* value 
measured on the CoRoT light curve with a very good preci- 
sion, and to get rid of the large uncertainty on the gravity 
estimate. It will be detailed in a next paper devoted to 
the fundamental parameters of the first CoRoT planet host 
stars, based on forthcoming UVES spectra. It is possible 
to assess that this star is quite an old main sequence star. 
The uncertainty on the mass, linked to the uncertainties on 
T e ff and [M /H] , will be significantly reduced using comple- 
mentary spectroscopic observations. We get M*=O.95M 
±0.15Mq. From the planet side, only the uncertainty on the 
stellar mass limits our knowled ge of the stellar and planet 
parameters (jBrown et al. 20011 ). 

7. Conclusion 

CoRoT-Exo-lb is a giant short-period planet orbiting at 
~ 5 stellar radii from a GOV star that seems to be metal- 
poor (further studies are necessary to confirm the metallic- 



Table 2. The star and the planet parameters. 



Star parameters 


CorotID 


0102890318 


RA [J2000] 


6 h 4g m l9 ir 


Dec [J2000] 


-3° 06' 07.78" 




13.6 


V [km/s] 


23.354 ±0.008 


K [m/s] 


188 ±11 


v sin i [km/s] 


5.2 ±1.0 


T eff [K] 


5950 ±150 


log g 


4.25 ± 0.30 


[M/H] 


-0.3 ±0.25 


Mi, [M ] 


0.95 ±0.15 


R* [R@] 


1.11 ±0.05 


Planet parameters 




k = Rp/Ri, 


0.1388 ±0.0021 


Tc [A] 


2454159.4532 ± 0.0001 


1 I'X , 

mI'^/r* 


0.887 ±0.014 


U+ 


0.71 ±0.16 


U- 


0.13 ±0.30 


P[d] 


1.5089557 ± 0.0000064 




4.92 ±0.08 


e 


(fixed) 


i [deg] 


85.1 ±0.5 


M P [M Jup ] 


1.03 ±0.12 


Rp [Riup ] 


1.49 ±0.08 


Pp [g cm" 3 ] 


0.38 ±0.05 


T P [K] 


1898 ±50 



ity). Its main characteristics are a large radius and a very 
low value of the mean densi ty that may be cons istent with 
a planet metal deficiency (B urrows et al. 20071 ). The low 
mean density could be interpreted as planet inflation unde r 
strong irradiation and high loss rate (jLammer et al. 2 003) . 
Similarly to HD209458b, the planet is too large to be re- 
produced by standard evolution models. Its radius could be 
explained by advocating an additional heat source (~ 1% 
of the incoming stel lar flux) or an inc rease of the 2000- 
3000K opacities (e.g. iGuillot et al. 20061) . The star and the 
planet are in strong tidal interaction with a Doodson con- 
stant of the same or der of magnitude as for OGLE-TR-56b 
([Patzold et al. 20"04l ). Due to its high photometric accuracy 
CoRoT will continue to help us improve our understanding 
of strongly irradiated planets. CoRoT-Exo-lb was detected 
early in the processing of the data. A lower level of noise 
with completely reduced data will open up the possibilities 
to search for outer planets, using transit timing variations, 
or to detect planet reflected light. 
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Table 1. Radial velocity measurements of CoRoT-Exo-lb 
obtained by SOPHIE associated to the system orbital 
phases (the large error bar of the second data point, due to 
acquisition problems, do not change significantly the orbital 
fitting). 



Julian 


Orbit 


RV 




date (d) 


Phase 


kms -1 


kms -1 


2454184.30498 


0.469 


23.2879 


0.0385 


2454185.30867 


0.134 


23.1284 


0.0600 


2454192.30397 


0.770 


23.5625 


0.0271 


2454197.32090 


0.095 


23.2862 


0.0309 


2454376.66501 


0.948 


23.5784 


0.0226 


2454378.66336 


0.273 


23.3730 


0.0324 


2454379.66472 


0.936 


23.6200 


0.0224 


2454379.67016 


0.602 


23.6907 


0.0229 


2454379.63151 


0.240 


23.3795 


0.0227 



